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Abstract 

Using pNRQCD and known results for the field strength correlator, we calculate the 
ultrasoft contribution to the QCD static energy of a quark-antiquark pair at short 



distances at N LO in a s . At the same order, this provides the logarithmic terms of 



the singlet static potential in pNRQCD and the log a s terms of the static energy. 



1 Introduction 



The ground state energy, E (r), of a static quark and a static antiquark sep- 
arated by a distance r is a key object for the understanding of the QCD 
dynamics. It is also a basic ingredient of the Schrodinger-like formulation of 
heavy quarkonium systems pQ. Its linear behavior at long distances is a signal 
for confinement [2], but also at short distances (r <C 1/Aqcd), where weak cou- 
pling calculations are reliable, it shows a non-trivial behavior. Indeed, when 
calculated in perturbation theory, infrared divergences are found starting at 
three loops [3]H] . These are due to the virtual emission of ultrasoft gluons with 
energy of the order E (r), which turn a color singlet quark-antiquark pair into 
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a color octet one and vice-versa. The proper treatment of the ultrasoft emis- 
sions requires the resummation of an infinite class of diagrams, which produces 
a non-analytic dependence on a s (typically logarithms of it). We shall focus 
here on this short distance behavior. 

The current knowledge of Eq{t) at short distance may be summarized as 
follows 
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where C F = T F (N* - 1)/N C , C A = N c , T F = 1/2, N c is the number of colors, 
A, = HC A /3 - 4T F n f /3 and ft = 34C|/3 - 20C A T F n f /3 - AC F T F n f are 
the first two coefficients of the beta function, rif is the number of (massless) 
flavors, 7^ is the Euler constant and a s is the strong coupling constant in the 
MS scheme. The one-loop coefficient ai is given by [5|6] 
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-T F nf, 
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and the two loop coefficient 02 by [7.8 jj9lfT0l 



/4343 2 7T 4 22 \ 2 A798 56 \ 

2= \ k i62 _ ~t y c(3) J ^ - br + t c(3) J 

- (f - 16C(3)) C F T F n f + (y Wf . (3) 

The logarithmic piece of the third-order correction was calculated in |llj . 
whereas the non- logarithmic piece a 3 has not been calculated yet. However, a 3 
is believed to be dominated by contributions which are known from renormal- 
ization group arguments [12]. If we write a 3 = 03 + 0^° , af G 3> a 3 . af G has a 
known expression in terms of the coefficients of the beta function and of those 
entering in the potential at lower orders (see [12], where a 3 = — 487r 3 x V3 , 
03 = 64cq). Estimates of 03 have been carried out using Pade approximations 
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[T2] and renormalon dominance [15)14)15(15] . which are consistent with the 
inequality above, and give similar results. The double logarithmic fourth-order 
correction may be obtained from [T7] , where higher-order contributions of the 
form a™ +3 log n a s were resummed using renormalization group techniques. It 
reads 

af = if!c5(-fc, + ir F% ). (4) 



The main result of this letter is the calculation of the logarithmic fourth-order 
correction to the singlet potential, a\\ we anticipate that it reads 



a\ = lQir 2 C 3 A 



_ /149 22, n 4 2 

+C A log 2 + -7T 2 

A V 27 9 6 9 



(5) 



The non-logarithmic piece remains unknown. 



A convenient method to calculate the logarithmic contributions to Eq. (CQ), 
which steam from the dynamics at the ultrasoft scale Eo(r), consists in inte- 
grating out from static QCD degrees of freedom at the soft energy scale 1/r and 
working within the effective field theory framework of pNRQCD |18.19](see 
[20] for a review). The quark-antiquark system may be in a color singlet or 
in a color octet configuration, which are encoded in color singlet, S, and color 
octet, O, fields in pNRQCD. At leading order in the multipole expansion, the 
integration of the soft energy scale gives rise to a singlet, V s (r;fi), and an 
octet, V (r; /i), static potential, which depend on r and a factorization scale \i. 
At next-to-leading order, two more "potentials" appear, V^(r; //) and V^(r; /i), 
which are the matching coefficients of the singlet-octet and octet-octet vertices 
respectively. At this order, the pNRQCD Lagrangian reads 



£pNRQCD — Aight 

+ J d 3 rTr|s t [i^-y,(r;/i)]S + Ot[iD -K(r;/i)]o| 

+V A (r; /i)Tr {O f r ■ #E S + S f r • #E O} 

+ ^ ;/i) Tr {O-tr -gEO + O^Or ■ ^e} 

+ ..-, (6) 



where £n g ht is the part of the Lagrangian involving gluons and light quarks, 
and coincides with the QCD one. The dots stand for higher-order terms in the 
multipole expansion. The static energy calculated from the above Lagrangian 
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has the form 

Eo(r) = V.(r; //) + 5 us (r, V s , V , V A , V B ,...;^), 



(7) 



where <5us( r , V s , V , Va, Vb, • • • ; A*) (^us f° r short) contains contributions from 
the ultrasoft gluons. V s (r; /x) and V D (r; fi) do not depend on \i up to N 2 LO [19] . 
The former coincides with Eq(t) at this order and the latter may be found in 
[2T] . The fact that the /i dependence of #us must cancel the one in V s (r; /j) 
is the key observation that leads to a drastic simplification in the calculation 
of the logo s terms in E (r). So, for instance, the logarithmic contribution 
at N 3 LO, which is part of the three-loop contributions to V s (r;fi), may be 
extracted from a one-loop calculation of <5us [TT1fT9] and the single logarithmic 
contribution at N 4 LO, which is part of the four-loop contributions to V s (r; fi), 
may be extracted from a two-loop calculation of 5usJZl 

In Sec. |2l we review the calculation of the third-order logarithmic term since it 
follows the same lines as that of the fourth-order one, which will be presented 
in Sec. El In Sec. 0], we conclude and discuss some applications of this work. 



2 Review of the third-order logarithmic correction 

In d dimensions, the order r 2 contribution due to ultrasoft gluons reads [TTirT9] 

2 2 00 

5 VS = -ijfTFVijZi I ^ e "^^ s) (°l Ea (O0(*,O)^'E 6 (O)|O). (8) 
c o 



<f)(y, x)^ 3 is the Wilson line in the adjoint representation connecting the points 
y and i by a straight line (t stands for (t, 0)). We will evaluate Eq. (jSJ) per- 
turbatively in a s . The dependence on a s , apart from the trivial g 2 factor, 
enters through (i) the V s and V a potentials, (ii) Va and (Hi) the field strength 
correlator of the chromoelectric fields. 

(i) The difference V — V s is given at leading-order by — L I , Note that 

2 r 

at leading and next-to-leading order V s and V Q only differ by an overall color 
factor. 

(ii) At tree level Va = 1. 



2 We denote N n LO, contributions to the potential of order a™ +l and N n LL, con- 
tributions of order a™ +2 log" --1 a s . 
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Fig. 1. Leading-order contribution to the field strength correlator. The double line 
represents the gluonic string, the circled cross and the springy line the chromoelec- 
tric field correlator. Note that our convention differs from the one in [23], where the 
gluonic string is either not represented, if no gluons emerge from it, or is represented 
by a dashed line. 

(Hi) The two-point field strength correlator 

2W*) = m{GiM<p{y^)^G b a*)} io> (9) 



can be parameterized in terms of two scalar functions T>(z 2 ) and T>i(z 2 ) ac- 
cording to 



+ (gn* Z v Z u ~ 9(JW Z v z \ ~ 9uX Z fi z ui + 9uui z fi z \) ' (10) 

where z = y — x [22]. In (jHJ), x and y only differ in the time component, hence 
z = t. Furthermore, in d dimensions, the chromoelectric component is given 
by 



(0\E a (y)cj>(y, x)^E 6 (x)|0> = V i0i0 (z) = -{d-l) 



V{z 2 )+V l {z 2 ) 
^dV^z 2 ) 



dz 2 



The leading-order contribution to the field strength correlator is given by the 
diagram shown in Fig. [TJ In d = 4 — 2e, the result is 



V?\z 2 )=^(N 2 -l) 



r(2 - e) 



TT 2 ~^-Z 2 ) 2 -^ 



V^{z 2 ) = 0. 



(12) 



Note that keeping e 7^ in the chromoelectric correlator provides a regular- 
ization for the integral over t in Eq. (JHJ). 

We now insert (i), (ii) and (Hi) into Eq. (J8|). Since at the ultrasoft scale 
t(V — V s ) ~ 1, the integral in t is performed without expanding the exponential 
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and making use of 



dtt n e- at = T{n ^ 1) . (13) 

a n+l V 1 







The final result reads: 



S vs = i^aj(l/r) f 4 - 2 log ^ + | - 2 log 2 ) , (14) 

24: r n \e /i 3 J 

1 1 

where — = 7e + log(47r). Note that the a s coming from the potential is 

evaluated at the soft scale 1/r, while the a s coming from the ultrasoft coupling 
is evaluated at the scale //. This will become relevant in the next section. The 
ultraviolet divergence in ( fl~4l) can be reabsorbed by a renormalization of the 
potential. In the MS scheme, in coordinate space, we have: 

K(r; //) - ZV s (r; »), Z = l + ^A^- a 2 s (l/r)-. (15) 



Since the static energy is \i independent, from the calculation above we infer 
that the logarithmic contribution to V s (r; //) at order a* must be 

5V s {r-n) = -CV^-^a 3 (l/r)logM, (16) 

1/ r 7T 



which added to the renormalized <5us contribution (from ( Tl4"j) ) gives the loga s 
term displayed in the third line of Eq. (pQ). This term was first calculated in 
|llj . where the cancellation between the IR cut-off of V s (r; fi) and the UV 
cut-off of the pNRQCD expression was checked explicitly by calculating the 
relevant Feynman diagrams in the Wilson loop. 

A comment is in order concerning the scheme dependence of the calculation of 
5us- This is not important if we are only interested in the logarithmic contri- 
bution, but it is if we wish eventually to combine our result with a (yet to be 
done) calculation of V s (r; jj) at N 3 LO and get the non- logarithmic pieces of the 
static energy right. We will assume that such a calculation will be done in mo- 
mentum space and that dimensional regularization and the MS scheme will be 
used to renormalize the UV divergences, like in the N 2 LO calculation [7118119] . 
The result will still be IR divergent when d — > 4, and the question is how one 
should proceed in order to combine that result with ours in a consistent wayf 3 ! 

3 Note that the MS subtraction of (|15|) in coordinate space is not equivalent to the 
MS subtraction in momentum space. 
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We propose to convert the (UV renormalized) momentum-space potential to 
coordinate space (in d dimensions) in that calculation, and together to use d- 
dimensional expressions for all the objects in our calculation, namely also for 
V s (r; /i) and V (r; fi) (Va(V; /i) remains the same in d dimensions). This guar- 
antees that the IR behavior of the regulated effective theory is exactly the 
same as the one of the fundamental theory. Had we expanded V Q — V s in Eq. 
([8]) we would have obtained zero, which means that the UV divergences, which 
remain after renormalization by the MS QCD counterterms (and by that of 
the color octet field wave function) in the effective theory, cancel exactly the 
IR divergences. Therefore, as a consequence of the fact that the IR behavior 
of the regulated effective theory is the same as the one of the fundamental 
one, the UV divergences in (TH|) cancel exactly the IR divergences in V s (r; /i), 
the /i dependence disappears, and the non-logarithmic pieces are correctly 
calculated. This procedure would be analogous to the one employed in [21] . 
Alternatively, one could use MS for the IR divergences of V s (r; //) in momen- 
tum space, work out the momentum space expressions for the (i-dimensional 
version of ffT4l) and make the MS UV subtraction accordingly. 

In the following section, we will use the same procedure employed here to 
obtain the next-to-leading IR logarithmic dependence of the static potential. 
That is the logarithmic contribution to the potential, which is part of the 
N 4 LO contribution. 



3 Fourth-order logarithmic correction 



Equation ([8]) does not rely on an expansion in a s , therefore it also provides 
NLO contributions to 5us- in fact, as we argue next, it provides the full con- 
tribution to this order. 

In principle, we may have diagrams with more insertions of the operators in 
([6]) and diagrams with operators of higher order in the multipole expansion 
that contribute to <5us at NLO. Concerning the former, for symmetry reasons 
we need at least two more operator insertions, which implies a suppression 
of «g with respect to the leading-order <5us- Concerning the latter, operators 
of higher order in the multipole expansion may be found in [25 26J. Their 
contributions are suppressed by al with respect to the leading-order <5us- To 
see this just recall that the ultrasoft fields (and derivatives acting on them) 

must be counted as Eq{t) ~ a s /r. Then, any insertion of the kind / dtr • E 

implies an a s suppression (with an extra a s suppression for any r'D-? acting 
on the chromoelectric field). For a given diagram, additional suppressions may 
appear due to the coupling constants in front of the chromoelectric fields. 
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Fig. 2. Static QCD diagrams for the leading-order matching of Va- The solid lines 
stand for a static quark and antiquark, the dashed line for a longitudinal gluon. 

The NLO contribution to <5tjs is then provided by Eq. (|Sj) evaluated at relative 
order a B . Since the dependence in a s enters through Va, V s , V and the chro- 
moelectric correlator, we need the 0(a s ) corrections to all these quantities. 
These will be given in the following two sections. Finally, in Sec. 13.31 we will 
obtain the fourth-order logarithmic correction to the potential. 



3.1 0(a B ) corrections o/Va, V s and V 



The 0(a s ) corrections to V s and V Q are well known. In particular, we have 

a„(l/r) 



V -V s = ^-a s (l/r) 
2 r 



l + (a x + 2 7e (3 ) 



An 



(17) 



The matching coefficient Va can be obtained by matching static QCD to pN- 
RQCD at order r in the multipole expansion. At leading order in a s , we have 
to calculate the diagrams shown in Fig. [2J They give the tree level result 
Va = 1. One may naively expect the first correction to be 0(a B ), but it is 
notH This becomes clear if we perform the calculation in dimensional regu- 
larization and in Coulomb eaue Indeed, the diagrams that we can draw 
at 0(a s ) correspond either to self-energy corrections or to iterations of the 
Coulomb potential, which are identical in the effective theory and hence do 
not contribute to the matching. Then, the first non-vanishing correction to 
the tree level result may possibly come from diagrams like the one in Fig. [31 
which is 0(a1) and, therefore, unimportant here. 



3.2 0(a s ) correction of the field strength correlator 



The 0(a s ) correction to the QCD field strength correlator was calculated in 
[23] . It is given by the diagrams in Fig. HI Here we need the expression in d 

4 The vanishing of the anomalous dimension of Va at one loop has been observed 
in [I7j. 

5 The potentials are independent on the gauge used in the matching. Therefore, we 
can use the most convenient one to do the computation. 



8 



I 



Fig. 3. Example of static QCD diagram that contributes to the next-to-next-to-lead- 
ing order matching of Va- 

dimensions because in (JSJ) the integral over t is singular. The <i-dimensional 
result for the a s correction is [27] 



V^(z 2 ) 



N C {N 2 C 



4tt 2 - 



7T 



2e 



r 2 (i 



z- 



2-2t 



n (N 2 - d Qs ^ ^ r 2 ri -e)(- 



2-2t 

■ 2 I 9i{t) 



(18) 
(19) 



with 



9i(e) 



3 + 8e - 6e 2 + 2e 3 B(-l + 2e, -2 + 2e) 



e(3-5e + 2e 2 ) 
6- 18e + 17e 2 -6e 3 
e 2 (3 - 5e + 2e 2 ) 
7>nj 1 - e 
N c e(3-2e)' 



3 - 2e 



(20) 



-2(l-e + e 2 ) 



£(-l + 2e, -2 + 2e) 
e(3 - 2e) 



(21) 



where 



B(u,v) = T(u)T(v)/r(u + v). 

Since the external points x and y are fixed, the divergences that we encounter 
in "Djoio coming from the expressions above should cancel against the vertex 
and gluon and octet field propagator counterterms. The counterterm for the 
vertex is zero, since, as seen in the previous section, the first correction to 
Va is of order a 2 . The counterterm for the gluon propagator is the usual 
one in QCD. The counterterm for the octet propagator coincides with the 
counterterm for the quark propagator in the Heavy Quark Effective Theory 
[28] but with the quark in the adjoint representation. We can represent the 
counterterm contributions by the diagrams of Fig. We have checked that: 
(i) the divergence coming from the first diagram in Fig. H]is canceled by the 
counterterm of the gluon propagator, (ii) the diagram (b) of Fig. HJdoes not 
give a divergent contribution (as one would expect from the fact that the 
gluons are attached to the external fixed points only) and (in) when we sum 
the remaining diagrams the divergence that we obtain is exactly canceled by 
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(a) (b) 




(g) (h) 

Fig. 4. Next-to-leading order contributions to the field strength correlator. The glu- 
onic string is represented by a double line. The shaded blob represents the insertion 
of the one-loop gluon self-energy. Symmetric graphs are understood for (c) and (d) . 
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Fig. 5. 0(a s ) counterterm diagrams for the chromoelectric correlator. The gluonic 
string (which comes from the octet propagator) is represented by a double line. 

the counterterm of the octet propagator. In the MS scheme, at 0(a s ), the 
contributions of the counterterms are given by 



V c - t -(z 2 ) = (22) 



ir 4vr 2 - e v '(-z 2 ) 2 ^e 



where in the brackets we have kept separated the —2 coming from the octet 

5 4 rif 

propagator counterterm from the h -Tp— ^- coming from the gluon propa- 

gator one. The renormalized d-dimensional result for the a s correction to the 
chromoelectric correlator is 



(i) 
mo 



-(3-2e) V {1 \z 2 ) + (-l + 2e)V^\z 2 ) 



+V c - t iz 2 ) + (-l + e)V c l t iz 2 ) 



(24) 



which, indeed, is finite for e — > 0. 



3.3 Calculation of the fourth-order logarithmic correction 



The results of the two preceding sections provide all the necessary ingredients 
to compute <5 US at NLO. Let us split Stj S as follows 

<5us = £(U>| 0(Qs) + Gv -v s \ 0{as) + ^i/r\ 0(as) , (25) 



where the first and second terms stand for the a s corrections to the field 
strength correlator and to the potentials respectively, and the last term ac- 
counts for the contribution induced by a change of scale in the N 3 LO calcula- 
tion. 
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First, we shall consider the contribution (1241) to the field strength correlator. 
After integration over t, which can be done using Eq. ({TBI , we obtain 



j(r 2 ) 

'{EE)\o( a 



(<*s) 



7T 



CI I 
r 



x 



A B „ , , 
^ + - + Ci log' 



C 2 log 



(26) 



with 



A- 
B- 
Cv 
C 2 -. 
D 



1 f 4T F n f 11 

21 " y° A 



-107>n/ + C A (6^ + 47 
1 / 4 TpUf 11 

el — s - + Y Ca 
i 



54 
1 
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4 T F n f (10-6 log 2) + C A ( -149 + 66 log 2-12 



7[- 



( 67 5 5 
TFn/ ( " y + q 7s + 5 log 2 - 2 log 2 2 - - log vr 



(27) 
(28) 
(29) 
(30) 



/ i 241 47 11 47 

+Ca { ~w ~ Y2 lE ~ 17 log 2 + y log2 2 + 12 log * - 12C(3) 



, 9 2 IE 2 2 , o , 71 1 

+ -7T - y7T -7T log 2 + y log 7T 



(31) 



Next, we display the contribution that we obtain if in (jSJ) we use the leading- 
order expression for the chromoelectric correlator but the 0(a s ) correction for 
V - V s : 



,(r 2 ) 

'Vo-V a \ 0{as) 



Cf—~ (ai + 27^0) 

7T 7T 16 r 



7T 

X 



1 , K-K 5 

— - log + - - log 2 

2e y ii I 6 



(32) 



The ultraviolet divergences in the expressions (1261) and f)32p come from the 
integration over time in (jSJ). They can be absorbed by a renormalization of 
the potential, analogous to ffT5l) . 

Finally, we obtain another contribution if in the renormalized version of (fill) 
we change a s (/i) to a s (l/r) (we want all a s evaluated at the scale 1/r): 
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,(r 2 ) 

M-+l/ r l0(a s ) 



a: 




) 



i „ 5 
+ log2- - 



(33) 



Adding up the renormalized versions of (|26|) and (1321) and equation (1331) . we 
obtain the contribution of 5us to E (r) at order a^. The complete calculation 
of Eq{t) at this order requires the knowledge of V s (r; fi) at the same order. 
However, to obtain the terms proportional to loga s it is enough to enforce 
E Q {r) to be independent of the factorization scale fi. This constrains the terms 
q,5 i g 2 r ^ anc j a 5 j Q g r ^ Q f s i n gi e t static potential to be 



Summing (I34p with (1251) provides the coefficients a^ 2 and a\ of the static 
energy E (r) given in Eqs. (jlj) and ([5]) respectively. 

Note that: (i) in order to cancel the \i dependence of the two double logarithms 
in <5us, l°g(r/i) l°g((V — K)/a*) an d l°g 2 ((l<^ — K)/aO> against the single double 
logarithm in 5V^, log 2 r/i, the coefficient of log(r/x) log((V^, — V s )/n) must be 
twice the one of log 2 ((V^, — V B )/fj). (ii) The coefficient of the double logarithm 
log 2 r\i in 5V S should coincide with the one obtained expanding the renorma- 
lization group improved static potential of The coefficients a\ 2 and 

a\ must be renormalization scheme independent] 6 1. We have explicitly checked 
that our result satisfies these requirements. 



4 Conclusions 

We have calculated the ultrasoft contribution to the QCD static energy of a 
quark- ant iquark pair at order oq . This is sufficient to obtain the logarithmic 
contribution to the pNRQCD singlet static potential at N 4 LO, which, in turn, 
provides the log 2 a s / r and log a s / r terms of the static energy of a quark- 
antiquark pair at distance r. The calculation heavily relies on effective field 



6 For a calculation of 5V S (r; /jl) in the subtraction scheme 1/e — 7^ + log7r we refer 
to [29]. 
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theory techniques and uses the result of Ref. [23j as a key ingredient. 

Possible applications of the result include precision comparisons with lattice 
data, heavy quarkonium spectra and t-t production near threshold. 

At short distances, the perturbative expression of the QCD static energy has 
been compared with lattice data at N 2 LO in [30T3T] and at N 2 LL in [14] . Our 
analysis provides a key ingredient for a N 3 LL analysis. 

Starting from the N 3 LO in a s , the quarkonium mass becomes sensitive to 
the ultrasoft scale, if the ultrasoft scale is assumed to be much larger than 
Aqcd [3"2"|33II24|I3"4] . Our result also provides an important ingredient for the 
calculation of the quarkonium mass at N 3 LL accuracy. 

Top-quark pair production near threshold, which will become an important 
production process at the ILC, is presently known at N 2 LO [32] ■ The cross 
section at N 2 LL (see e.g. [351137] ) and at N 3 LO (see e.g. [MBH]) is computed 
presently by several different groups. Our result will contribute to the cross 
section at N 3 LL. The third-order renormalization group improved expression 
will be needed to resum logarithms potentially as large as the N 3 LO and reduce 
the scale dependence of the cross section. 
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